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Abstract 

The production of fast neutrons (1 MeV - 1 GeV) in high energy muon-nucleus 
interactions is poorly understood, yet it is fundamental to the understanding of the 
background in many underground experiments. The aim of the present experiment 
(CERN NA55) was to measure spallation neutrons produced by 190 GeV/c muons 
scattering on carbon, copper and lead targets. We have investigated the energy 
spectrum and angular distribution of spallation neutrons, and we report the result 
of our measurement of the neutron production differential cross section. 
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1 Introduction 



Neutrons are an important source of background in many low rate under- 
ground experiments. Detectors such as the Palo Verde detector [1], the 
Kamiokande detector [2], the CDMS detector [3], the Gallex detector [4], and 
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most of the other underground experiments are concerned with neutron back- 
ground. It is therefore important to understand the different neutron produc- 
tion mechanisms in the environment of such experiments. 

There are two sources of neutrons in underground laboratories. The first one 
is the natural radioactivity such as {a,n) reactions and spontaneous fission, 
due to uranium and thorium traces in the rock or in the materials used in the 
experiments themselves. The second source is the interaction of atmospheric 
muons, produced by cosmic rays, with matter, such as the surrounding rock, 
the different shieldings or the detector itself. Propagating through matter, 
muons lose energy mainly continuously by ionization. In addition, they lose 
energy in discrete bursts along their trajectory by bremsstrahlung, direct pair 
production and nuclear interactions [5]. The goal of the CERN NA55 exper- 
iment was to study neutron production in muon-nucleus interactions, which 
at present is poorly understood. In this paper we report the result of our 
measurement of the cross section for fast neutron production by interaction of 
high energy muons off carbon, copper and lead nuclei. Specifically, we inves- 
tigated the energy spectrum and angular distribution of spallation neutrons 
produced by the inelastic scattering of 190 GeV/c muons on each target, via 
the time-of-flight method. 



2 Interaction of high energy muons with matter 

The total neutron yield measured at different depths has been reported in 
references [1] and [6,7]. However, very little is known about the energy and an- 
gular distributions of the neutrons, even though the neutron energy spectrum 
in a scintillator has been determined in [7] for neutrons which penetrated at 
least 1 m from the muon track. Moreover, one expects that many neutrons are 
produced in secondary processes following the primary muon-nucleus interac- 
tion and it is difficult to separate the two effects in the inclusive measurements 
performed underground. 

To describe the process theoretically, one assumes that the muon scattering oc- 
curs primarily at very small angles. Such a process can be viewed as produced 
by a beam of nearly real equivalent photons [8]. Even with this drastic simpli- 
fication, it is difficult to obtain reliable estimates for actual yields and spectra 
of the neutrons. The difficulty lies in the estimation of the equivalent photon 
flux, as well as in the shortage of experimental nuclear photo-disintegration 
data for high energy photons. 
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3 Experimental appeiratus 

3.1 Experimental set-up 

The layout of our NA55 experiment is similar to a CERN LEAR experiment [9] 
in which neutrons and pions were measured following anti-protons stopping 
in the target. Another group, E665 at Fermi-lab [10], measured neutrons pro- 
duced by deep-inelastic scattering of 470 GeV/c muons. However, the latter 
experiment measured neutrons in a small energy range from 1 to 10 MeV, with 
high-energy deep inelastically scattered muons from nuclear targets. In the 
frame of the CERN NA54 experiment, measurements were done on a carbon 
target to study the production of radioactive isotopes in scintillation detectors 
by 100 and 190 GeV muons and their secondary shower particles [11]. 

The M2 muon beam at the CERN-SPS was well suited for our investigation. 
At 190 GeV/c, the energy of the SPS muons is similar to the mean energy of 
cosmic-ray muons at many underground experimental sites. The beam trans- 
verse size is 2.2 cm (FMWH). Our experiment was performed with a cylindrical 
target of 75 cm length and 8 cm diameter for the graphite, 25 cm length and 10 
cm diameter for the copper, 10 cm length and 20 cm diameter for the lead tar- 
get. Three neutron detectors (see section 3.2) were placed at 45 (Nl), 90 (N2) 
and 135 (N3) degrees relative to the muon beam, for a target-detector distance 
of 2.20 m (sec Fig. 1). A beam counter H6 (from NA47 experiment [12]) con- 
sisting of plastic scintillators is placed directly in the muon beam and is thus 
exposed to rates of about 20 MHz, during the two-seconds beam-on phase. To 
reduce the detector rate to a manageable level, the beam counter is segmented 
into 4 rings, each divided in 16 counters thus reducing the rate per phototube 
to at most 1 MHz in the center. Several large plastic scintillators are placed 
upstream from the neutron detectors and are used to veto muons originating 
from the beam halo. An iron hadron absorber was placed far behind the target 
along the beam axis. 



3.2 The neutron detectors and electronics 

Each detector is a 20 x 20 cm cylindrical glass vessel, equipped with two pho- 
tomultipliers (Phihps XP4512 photomultipliers with fast 1.1 nsec RMS pulse) 
and filled with Bicron 501 A liquid scintillator which possesses excellent pulse 
shape discrimination properties for neutron identification [13]. Thin plastic 
scintillators (7 mm) are positioned in front of and around each neutron detec- 
tor to tag charged particles, such as protons and pions. They are also used to 
veto the cosmic muons. 
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The particle energy is measured via the time-of-flight (TOF) from its produc- 
tion site in the target. The signal of the neutron detector provides the TOF 
start, and the signal of the segmented beam counter, located 6 m upstream 
from the target, is delayed to provide the TOF stop. The analog signals of 
each of the neutron counters N are coded using 3 different CAMAC ADC's 
digitizing respectively the full analog pulse, the pulse rise and the pulse tail. 
The pulse itself provides the three ADC gates by means of a low-threshold 
discriminator, which also provides a CAMAC multi-hit TDC with the TOF 
start signal. The analog signals of each of the plastic scintillators S are coded 
using 15 different CAMAC ADC's. 



4 Neutron detector efficiency 



The neutron detection efficiency was simulated using the GCALOR Monte 
Carlo package [14] . Neutrons were generated with energies from 2 to 1000 MeV. 

The efficiency for neutrons entering the detector was obtained as a function of 
both the detector energy threshold and the neutron energy. The geometrical 
acceptance of the detector, the quenching effect and each detector threshold 
were taken into account. The calculation of scintillator quenching was taken 
from Ref. [15], as the amount of light depends on the particle type and is 
not directly proportional to the deposited energy. Above a neutron energy of 
50 MeV, the electron equivalent energy in the scintillator is assumed to be a 
linear function of the proton energy [16], with a proportionality constant of 
0.83. Below 50 MeV, we apphed the following formula: 



Te = 0.83 Tp - 2.82 



1.0 - e ( - (1) 



where Tg was the electron kinetic energy in MeV and Tp the proton kinetic 
energy in MeV. 

As will be shown later, the neutron energy spectrum in each counter is peaked 
at low energies. Therefore, the detection efficiency is strongly dependent upon 
the neutron detectors thresholds. The thresholds are determined from the 
energy loss spectrum of cosmic muons in the scintillator cells. Vertical muons 
(cosmics) are selected by requiring a coincidence signal above pedestals in 
the top and bottom scintillation counters during beam-off runs. Because of 
the hardware thresholds, these spectra show a clear cut-off at low energy. 
The neutron detector thresholds are then extracted by comparing a simulated 
energy spectrum of cosmic muons with the measured one. The Monte Carlo 
spectrum of cosmic muons N{E) is then convoluted with a gaussian response 
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function, R{E, H) for the detectors: 

E 

N{H) = j R{H, E)N{E)dE (2) 



where N(H) is the measured spectrum in ADC counts H. The response func- 
tion parameters arc of the same form as the ones used by Arneodo et al. [13] 
who have described and tested the response of similar Bicron 501 A hquid 
scintillator cells. Applying this procedure, we obtain thresholds in electron 
equivalent energy of 4.8 ± 0.4 MeV for Nl, 5.5 ± 0.3 MeV for N2 and 6.9 ± 0.3 
MeV for N3. In neutron or proton recoil energy they are 8.7±0.7 MeV, 9.6±0.5 
MeV and 11.4±0.5 MeV for Nl, N2 and N3 respectively. Fig. 2 shows a com- 
parison of the measured and simulated, N{H), energy spectra of cosmic muons 
for the N3 detector. The errors on the thresholds are assigned to be 2 channels 
(see Fig. 2) and are overestimated. 

Finally, Fig. 3 shows the Monte Carlo simulation of the neutron detection 
efficiency as a function of both neutron energy and energy threshold. Neu- 
trons were assumed to be detected only if recoil protons ionize heavily and 
stop in the scintillator producing delayed light and thus making pulse shape 
discrimination possible. 



5 Data analysis 

The results of the analysis of our data from the different targets is presented 
below. 



5. 1 Background 

The challenge of this experiment is to discriminate neutrons from the abundant 
bremsstrahlung photons. Furthermore, four other background sources may 
contribute to the measured spectra: charged particles from the target, cos- 
mic muons, muons from the beam and ambient neutrons. Cosmic muons and 
charged particles can be rejected using the plastic scintillators. Muons from 
the beam halo are ehminated by the dedicated veto counter (HV, see Fig. 1). 
Surrounding neutrons cannot be eliminated, but their contribution can be es- 
timated from dedicated empty-target runs. Finally, neutrons and photons are 
identified by a pulse shape discrimination (PSD) of their signal in the liquid 
scintillator. 
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5.2 Pulse shape discrimination 

The time-development of scintillation light depends on the nature of the ion- 
ising particle. In scintillators such as Bicron 501 A higher ionisation density 

due to slow charged particles, such as recoil protons from neutron scattering, 
leads to the excitation of longer lived states and thus to a slower component 
in the light pulse. Measuring the signal amplitude at 2 different time inter- 
vals allows to visuahze this in a PSD plot [17]. The complete PSD method 
is explained in Ref. [18]. Fig. 4 shows the amplitude obtained in a delayed 
window as a function of the amplitude measured in a prompt one (PSDl). 
Based on the same principle. Fig. 5 is obtained by comparing the amplitude 
in the prompt window with the total signal amplitude (PSD2). In these two 
different PSD representations, we can distinguish a 7-like curve containing 
bremsstrahlung photons and cosmic and remaining beam halo muons, and a 
neutron-like zone containing all the neutrons and the charged hadrons from 
the target. The 7-likc events are easily selected from the time-of-flight since 
they are prompt with respect to the beam muons. Then we calculate for each 
event the real distance (in units of channel number) between the correspond- 
ing amphtudes and the fitted 7-like curve for PSDl (distance 1) and PSD2 
(distance 2). We can then represent distance 1 as a function of distance 2, for 
all events including the charged particles. Figure 6 shows the events left after 
the cuts which select only neutral particles. A clear separation can be seen 
between photons and neutrons. A cut is applied to select neutron like events 
achieving nearly complete suppression of the 7-like events. The acceptance is 
calculated by applying a gaussian fit on the projected neutron and gamma 
peaks, and evaluating the relative area under the gamma peak above thresh- 
old. The acceptance was calculated separately for each detector and target. 
This was necessary because of the hmited stabihty of the PSD, and because of 
the differences in the relative strengths of the neutron and gamma peaks from 
case to case. The global cut acceptances, reflecting these variations between 
targets and detectors, are shown in Table 1. The accepted neutron events are 
then used for the time-of-flight analysis. 



5.3 Time-of-flight analysis 

Before being converted into energy distributions, the selected neutron time- 
of-flight spectra must be corrected for off-target events and for random TOF 
stop signals. The first background arises from triggers correlated with beam 
particles, as neutrons bouncing on the ground or walls, for example. Their 
contribution can be deduced from 'empty target' runs properly normalized 
to the incident flux. The second background component is due to multiple 
TOF stop signals given by randomly incoming muons hitting the H6 detector 
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scintillators. As their time distribution is flat (see Fig. 7), they can be easily 
subtracted. 

The calibration of the time spectrum is obtained from the gamma peak. The 
time resolution slightly depends on the target. As an example, it is 1.0 ns, 4.6 
ns and 3.8 ns for Nl, N2 and N3 respectively for the carbon target. In this 
case we measure maximum neutron energies of around 1 GeV, 250 MeV and 
310 MeV for Nl, N2 and N3 respectively. The time resolution determines the 
energy resolution, which is shown as function of energy in Fig. 8. The energy 
distribution is computed from the TOF distribution, and corrected for the 
detector efficiency as described in section 4. 



6 Results 

6.1 Additional acceptance corrections 

The charged particle rejection using the scintillators S leads to a loss of good 
neutron events. The corresponding acceptances are determined from the com- 
parison of neutron samples with and without scintillator cut for each neutron 
detector, and strongly depend on the scintillator efficiencies and thresholds. 
In addition, the 85 Hz trigger rate for all detectors above a threshold of 50 
MeV leads to a dead time of 2.1%. All the efficiencies taken into account in the 
calculation of the neutron production cross section are summarized in Table 1. 

6.2 Production cross section and energy spectrum 

The neutron production cross section was calculated using the following for- 
mula: 



where dn is the number of selected neutrons corrected for all acceptances. All 
parameters are given in Table 2. Figures 8, 9 and 10 show the differential neu- 
tron production cross section as a function of the neutron energy, including all 
efficiency corrections, for the three angles, and for the different targets. The 
energy integrated values from the corresponding threshold up. are gathered 
in Table 3. It is worth noting that the neutron energy spectrum is displaced 
towards low energies when going from forward to backward angles. As a con- 
sequence, no neutrons are detected above 70 MeV at 135°. Our measurements 




(3) 
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show that the angular distribution is forward peaked, with the exception of 
the lowest energy bin. This effect can be understood as due to the recoil of the 
source nucleus combined with an isotropic distribution of the emitted neutron 
at the lowest energies. One also observes that the angular distribution gets 
wider with increasing atomic number. Generally our angular and energy dis- 
tributions agree with the FLUKA simulation performed by Wang et al. [19]. 
The slope of the neutron energy spectrum for carbon is in a qualitative agree- 
ment with the slope ~ (integrated over all angles) found in Ref. [7]. 

Figure 1 1 shows the integrated cross section as a function of the atomic number 
for the three different angles and three different targets. We note that the cross 
section increases significantly with the atomic number. 



7 Conclusion 



We investigated the production of spallation neutrons obtained from 190 GeV/c 
muons scattering on graphite, copper and lead targets. The neutrons were ob- 
served by liquid scintillator detectors, allowing background rejection by means 
of pulse shape discrimination. The neutron energy distribution was determined 
via time-of-flight. The 190 GeV muon energy corresponds to the mean energy 
of cosmic-ray muons at underground experimental sites of about 2000 meters 
water equivalent depth. 

It should be noted, however, that in the present experiment only neutrons 
associated with the primary muon-nuclear spallation process are detected, 
while in the underground detectors the total neutron yields are measured. Thus 
the present experiment is not dependent on the subsequent neutron transport 
and multiplication. In view of this, and also in the view of the difference 
in the muon spectrum (narrow vs. broad distribution) direct comparison is 
difficult. Nevertheless, the cross sections quoted in Table 3 are in a qualitative 
agreement with the yields underground and, as mentioned above, their shape 
agrees with the FLUKA simulation. 

Results on neutron angular distribution and energy distribution were obtained 
for the first time. The differential cross section as a function of the neutron 
energy were obtained for 45°, 90° and 135° production angles with around 
15% accuracy. 
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HV 




Fig. 1. NA55 detector layout. N: Bicron 501A liquid scintillator detectors with 
Pulse Shape Discrimination. S: thin plastic scintillator counters for charged particle 
identification, in front and around each N detector. H6: NA47 segmented high rate 
beam counter. HV: beam muon halo veto, thin plastic scintillator. Arbitrary figure 
scale. 



Table 1 



Acceptance of the cuts 



Target 


Analysis cuts 


Nl (45°) 


N2 (90°) 


N3 (135°) 


Carbon 


PSD cuts 


87% 


66% 


64% 




Plastic scintillators 


67 % 


82 % 


66 % 


Copper 


PSD cuts 


48% 


91% 


26% 




Plastic scintillators 


69 % 


85 % 


88 % 


Lead 


PSD cuts 


50% 


28% 


28% 




Plastic scintillators 


57 % 


84 % 


87 % 




Life time 




97.9% 
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Table 2 



Summary of calculated cross-section parameters 









Angular acceptance dfl 




(5.2 ± 0.6) X 10-^ sr 


Muons on target (/) 


C 


1.63 X 10^^ muons 




Cu 


1.84 X 10^^ muons 




Pb 


1.93 X 10^^ muons 


Target density N=pL 


C 


(8.52 ± 0.11) X 10^4 atoms/cm^ 




Cu 


(2.12 ± 0.03) X 10^4 atoms/cm2 




Pb 


(3.30 ± 0.04) X 10^^ atoms/cm^ 


Total number of neutrons, 




45° 100500 neutrons 


above threshold, corrected 


C 


90° 53510 neutrons 


for the detector efficiency 




135° 29395 neutrons 



45° 277960 neutrons 
Cu 90° 135545 neutrons 
135° 148955 neutrons 



45° 251360 neutrons 
Pb 90° 316820 neutrons 
135° 105060 neutrons 



Tabic 3 



Dilfcrciilial iic}ulroii production cross-s(;cli(jii 



Differential cross section 


target 


barn/sr 




neut//x-sr-g-cm ^ 




^(45° ±1°) 

^(90° ±1°) 
^(135° ±1°) 


C 


(0.21 ± 0.04) X 10" 
(0.09 ± 0.02) X 10" 
(0.06 ± 0.01) X 10" 


-3 

-3 

-3 


(1.06 ± 0.20) X 10" 

(0.46 ± 0.08) X 10" 
(0.31 ± 0.06) X 10" 


-5 

-5 

^5 


^(45° ±1°) 
^(90° ±1°) 
^(135° ±1°) 


Cu 


(2.00 ± 0.40) X 10- 
(0.80 ± 0.15) X 10" 
(0.83 ± 0.15) X 10- 


-3 
-3 
-3 


(1.90 ± 0.38) X 10" 
(0.76 ± 0.14) X 10- 
(0.79 ± 0.14) X 10- 


-5 
-5 
-5 


^(45° ±1°) 
^(90° ±1°) 
^(135° ±1°) 


Pb 


(13.53 ± 3.00) X 10' 
(11.67 ± 2.38) X 10" 
(3.75 ± 0.72) X 10" 


-3 
-3 

-3 


(3.93 ± 0.87) X 10- 
(3.39 ± 0.69) X 10- 
(1.09 ± 0.21) X 10- 


-5 
-5 
-5 
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100 200 300 400 500 6C 

H (ADC counts) 



Fig. 2. Measured (error bars) and simulated (histogram) energy spectra of cosmic 
muons for the N3 detector. The lower energy threshold, on the left, is clearly visible 
and well determined within 2 channels and corresponds to 6.9 it 0.3 MeV electron 
equivalent energy. 
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Fig. 3. Simulation of the neutron detection efficiency as a function of neutron energy 
and thrcsliold (in electron equivalent energy), taking into account the geometrical 
acceptance of the detector and the quenching effect, for the angles 45° (Nl), 90° 
(N2) and 135° (N3). Polynomial interpolations (dotted, solid and dashed lines) have 
been obtained from the simulated efficiencies (open circles, solid circles and open 
squares) and used for the efficiency corrections. The bump at 700 MeV is the result 
of an approximation. 
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Fig. 4. First pulse shape discrimination 
(PSDl) for the Nl detector and the car- 
bon target. The signal is measured in two 
different time intervals : prompt signal 
amplitude vs. delayed signal amplitude. 



Fig. 5. Second pulse shape discrimina- 
tion (PSD2) for the Nl detector and the 
carbon target. Prompt signal amplitude 
vs. total signal amplitude (different bin- 
ning) . 



gammas 




neutrons 



■300 - 600 -400 -200 200 400 600 800 1000 1200 



Distance 2 (channel number) 



Fig. 6. Distance 1 vs. distance 2 for the Nl detector and the carbon target. These 
distances were calculated between every point and the 7-like events for PSDl and 
PSD2. Only neutral particles are selected. 
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Fig. 7. TDC spectrum corresponding to the neutron distribution dn/dt for the 
carbon target. The undershoot for times immediately before the 'prompt' signal is 
caused by the removal of the gamma peak. The background due to uncorrected 
beam muons is flat. 
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3 



Neutron energy (MeV) 



Fig. 8. Neutron energy spectrum for the carbon target, including all the corrections. 
The spacing between two ticker marks on the upper horizontal line gives the energy 
resolution (FWHM) at the corresponding energy, derived from the time resolution, 
for the detector Nl. 
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Neutron energy (MeV) 



Fig. 9. Neutron energy spectrum for the copper target, including all the corrections. 
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Fig. 10. Neutron energy spectrum for the lead target, including all the corrections. 
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Fig. 11. Integrated cross section vs atomic number for the three different angles. 
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